Ground state of ID bosons with delta interaction: 
link to the BCS model 

M.T. Batchelort§, X.W. Guant and J.B. McGuireJ 

t Department of Theoretical Physics, RSPSE, and 

Centre for Mathematics and its Apphcations, MSI 

The Austrahan National University, Canberra ACT 0200, Australia 

I Physics Department, Florida Atlantic University, Boca Raton, Florida 33431 

Abstract. The Bethe roots describing the ground state energy of the integrable ID 
model of interacting bosons with weakly repulsive two-body delta interactions are seen 
to satisfy the set of Richardson equations appearing in the strong coupling limit of an 
integrable BCS pairing model. The BCS model describes boson-boson interactions 
with zero centre of mass momentum of pairs. It follows that the Bethe roots of 
the weakly interacting boson model are given by the zeros of Laguerre polynomials. 
The ground state energy and the lowest excitation are obtained explicitly via the 
Bethe roots. A direct link has thus been established, in the context of integrable ID 
models, between bosons interacting via weakly repulsive two-body delta-interactions 
and strongly interacting Cooper pairs of bosons. 
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There has been a revival of interest in the exactly solved ID model of interacting 
bosons pp. The conventional description of the model is N interacting bosons governed 
by the Hamiltonan 

^ = -E9l2+2c 6{x,-x,) (1) 

i=i * i<«<i<A^ 

and constrained by periodic boundary conditions to a line of length L. Here 2c is the 
strength of the two-body delta interaction, with h = 2m = 1. The wave functions are 
given in terms of the Bethe ansatz by 

^{xi, ■ ■ ■ , xtv) = ^ A{p) exp(i ^ kp^Xj) (2) 
p j 

in the region xi < X2 < x^ < ■ ■ ■ < xn- The summation extends over all permutations 
p of momenta {kj} and A[p) are coefficients depending on p. The eigenvalues are given 
by £^ = J2f=i ^"j where {kj} satisfy the Bethe equations 

exp(iA:,L) = - n ^^^1^ for j = l,...,iV. (3) 

■^■^ k,; — kf ~ IC 

The Bethe roots kj are known to be real for repulsive interactions, c > 0, for which the 
ground state energy and excitation spectrum have been studied extensively. || 

The revival of interest in this model has been inspired by recent experimental 
and theoretical work on low-dimensional trapped boson gases at ultracold temperatures 
[HI 13 Uni UH 1121 Uni UH USI; particularly with regard to highly elongated traps which 
give possible realizations of a ID quantum gas as well as Bose-Einstein condensates. In 
a remarkable recent experiment, a ID quantum gas was created in an optical lattice 
in which the interactions between trapped ultracold Rb atoms were modified to bring 
about a continuous passage from the weakly interacting regime to the strongly correlated 
Tonks-Girardeau gas (Ref. jH] and refs therein). The key parameter is essentially the 
ratio 7 of interaction to kinetic energy. In this way the pronounced fermionic behaviour 
of the Tonks-Girardeau gas was observed in the strong coupling regime jH]. 

Here we consider asymptotic solutions to the Bethe equations of the interacting 
boson model in the strong and weak coupling limits. We see that the momentum 
density distribution agrees with the observed profiles of the ID interacting quantum 
gas [Hj. Moreover, we find a link between the Bethe equations (jHI) in the weak coupling 
regime and the Richardson equations ^H] for the BCS boson pairing model in the strong 
coupling regime I18| IT^. The BCS model describes boson-boson interactions with 
zero centre of mass momentum of pairs. In this way the Bethe roots for the ground 
state energy of the weakly interacting Bose gas are characterized by either the roots of 
Hermite polynomials or appropriate Laguerre polynomials. In the other direction, we 
thus see an explicit connection between strongly interacting Cooper pairs described by 
the BCS boson model and bosons interacting through weakly repulsive two-body delta 
interactions. We begin by looking at the strong coupling limit. 

II See, e.g., Refs PEIEIIIEIEIIZI and references therein. 
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Strong coupling limit 

In the Lc = oo limit, for the dilute gas {L > N) there is a well known connection to 
non- interacting fermions, known as the Tonks-Giradeau gas |2D]- The Bethe equations 
(121) reduce to exp{ikjL) = (—1)^"^, with solutions 2£'7i/L with £ integer for odd N and 
half-odd-integer for even A^. Corrections to order 1/c can be readily obtained from 
the asymptotic solutions of the Bethe equations Q- Define A = Lc/N. The roots are 
real and symmetric about the origin, with |±/c2m-i; = 1? ■ ■ ■ ? y} for even A^, where 
ki = X (-'■ ~ f)- ^'^^ '^^^ roots are zero and [^±.k2m-,fn = 1, ■ ■ ■ , "^^^j- The 

ground state energy follows directly from these asymptotic solutions, with 

where p = N/L. This result coincides with that of the perturbation theory approach 
(see, e.g. Ref [21j). At zero temperature, the chemical potential is given by /i = 
vr^p^ — 167r^p^/3c. The two-body correlation function g2 = An'^pl/'y'^ follows from the 
ground state energy and the Hellmann-Feynman theorem, with g2 giving the rates of 
two-body inelastic processes in ID trapped boson gases ^3]- Moreover, we see that 
the momentum density distribution for the ID A^-body interacting bosons is flat due 
to the symmetric equal-spacing distribution in momentum space. This qualitatively 
coincides with the experimental momentum profile of the Tonks- Girardeau gas in the 
strong coupling limit, for example, 7 ~ 204.5 [Hj. In this regime, the strong repulsive 
delta potential prevents the bosons from occupying the same position. Therefore, they 
are spread out and extend to a larger region in momentum space than in the case of 
weakly interacting bosons. 

Weak coupling limit 

The situation in the weak coupling limit is far more interesting. The ground state 
is an analytic function for c > 0. Considering the limit Lc << 1 and after some tedious 
case by case calculations for A^ = 2, 3, . . . , 12, we find that the Bethe roots satisfy the 
set 

of nonlinear algebraic equations. Here the summation excludes j = i and dj = 
0, ±1, ±2, . . . denotes excited states for fixed A^. The ground state has zero total 
momentum, with dj = for j = 1, . . . ,N. These equations are closely related to 
Stieltjes problems [221 • The ground state energy per particle 

So ^ cjN-l) 
N L ' 

follows directly from dj}. The nonlinear algebraic equations (0) for the ground state 
were found by Gaudin [21 , who showed that the kj are roots of Hermite polynomials of 
degree A^, namely Hjq{k) = 0. 
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We are also interested in the lowest excited state, which has total momentum 
k = Without loss of generality, we may take di = 1, d2 = ■ ■ ■ = d^ = 0. The lowest 
excitation energy per particle 

Si _ c{N - 1) , 2n (2tx ^ 2c{N - I) ^ ^{N -2f \ 



N L LN \L L L27r ^ 

follows from solving equations (jS)). Correspondingly, the largest Bethe root is fci = ^ + 
2c(jv-i) _|_ clearly seen that the energy gap vanishes in the thermodynamic 

limit. Some numerical exploration reveals that the excitation energy per particle (|7j) is 
very accurate for arbitrary number of bosons. Other excitation energies depend on the 
assignments of dj, with total momentum k = "Y^^ 2djT^ / L. 

The connection between the Bethe roots for the ground state energy and the roots 
of Hj^{k) provides a systematic way for studying quantities such as correlations and 
the momentum distribution function. The normalized momentum density distribution 
is given by the semi-circle law jSj 

nik) = ^^{Acp-e)K (8) 

We show the momentum density distribution in FigureHobtained by numerically solving 
the Hermite polynomials of degree = 60 and 100 for different values of c. The results 
fit very well the analytical expression (jSI). A similar semi-circular momentum density 
distribution has been found in the ferromagnetic Heisenberg chain j2S] • We see that the 
distribution at small momentum is rather flat, with an almost linearly decreasing region 
at large momentum. The stronger the interaction strength, the larger the momentum 
distribution region. This reveals a significant signature of the ID boson gas in the 
weakly repulsive limit Lc « 1. Remarkably, this behaviour was recently observed in 
the experiment for weakly interacting bosons with 7 ^ 0.5 fS]. If c = all the particles 
condense in the ground state at zero temperture. 

In order to link the ID boson gas with two-body delta interactions to the BCS type 
of pairing models, we consider first an even number of bosons, with A^ = 2M. We also 
find, to first order in c, that the Bethe roots for the ground state are of the form 



ki,2 = ±V El, /c3^4 = ± V^a, . . . , k2M-i,2M = ±V Em, (9) 
where the Ei satisfy the equations 

_L ^' 2 1 

i=i ^ 

for i = 1, . . . , M, where the summation excludes j = i. 

Now similar equations have arisen in a number of contexts jT3 I22|- Of particular 
interest here is the connection between Eq. (fTUI) and Richardson's equations for the 
BCS pairing model fl^ ES] in the strong coupling limit ^^I- However, the precise 
link between the integrable boson model with weakly repulsive delta interaction and 
the standard BCS model with strong attractive pairing interaction is quite subtle. In 
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Figure 1. The density n(fc) of momentum distribution (normalized) for TV = 60 
and 100 particles: the density distribution function sensitively depends on interaction 
strength, particle number and the length L for the ID boson gas with repulsive delta 
interaction in the limit Lc << 1. 



terms of the Cooper pair energies Ei participating in the scattering process, these latter 
equations are 

1 ^ y^' 2 _n' ^^^^ 

^Bcsd ^ Ei — Ej Ei ' 

for i = 1, . . . , r. Here n' = n — 2{m — r), where n is the number of unblocked energy 
levels, m is the number of Cooper pairs and r is the number of nonvanishing roots to 
order Abcs; where Abcs is the superconducting coupling constant, with d the mean level 
spacing. 

In the BCS model, the parameters {Ei} in Eq. (fTT|) represent the energies of r 
Cooper pairs participating in the scattering. Now the ground state of the BCS model 
in the strong coupling limit has r = m {n' = n), the first degenerate group of excited 
states corresponds to r = m — 1, etc |17,. Direct comparison between Eqs. (fTUj) and 
|TT|l gives r = M and n' = — |, leading to the ground state where total Cooper pairs is 
m = r = M. However, in the strong coupling limit, all of the energy levels collapse into 
one multiply-degenerate level, namely n = 1. 

Obviously, this state appears to lie outside the physical regime of the BCS model 
[T7] due to the presence of the negative fractional quantum number n' = — |. The main 
reason for this is that in the standard BCS pairing interaction, the electrons are paired 
into singlets rather than into triplets with zero centre of mass momentum of pairs. 
Therefore the pairing interaction term contains only a pair of attractive electrons with 
opposite spin and momenta so that the degeneracy at each energy state is a doublet, 
with level degeneracy Q = 2 j23]. In addition, the commutation relation between BCS 
pairs (hard core bosons) is not the same as that for bosons. 
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Before addressing the correspondence, we can nevertheless make use of the 
observation fTT] that the nonvanishing roots of the BCS Eq. (fTT|) are given by the 
zeros of an associated Laguerre polynomial. It follows that the Bethe roots of the 
weakly interacting boson gas in Eq. ()10|) are given in terms of the associated Laguerre 
polynomial L^{x) by 

LJ = 0, (12) 



. 2c 

for i = 1 M, where 



^"^"^^ ^^Jn-rl)!(A; + m)!m!^"'- ^^^^ 

Explicitly, for M = 1 (A^ = 2), = c/L and for M = 2 (A^ = 4), = 
c(3 - y6)/L, E2 = c(3 + VQ)/L. 

For an odd number of bosons, = 2M + 1, we find the ground state Bethe roots 
to be given by 

h = 0, ^2,3 = ±a/^, h,5 = ±\/E^, • • • , k2M,2M+l = ±\/E^I, (14) 

where the Ei satisfy the equations 

A^' 2 L 3 

^ Ei-E, ~Yc^~2Ei' ^^^^ 
j=i J 

for i = 1, . . . , M. Comparing again with Eq. n' = — | with again r = m = M. In 
this case the Bethe roots are given by 

4 = 0, (16) 



V 2c 

for i = 1, . . . , M. The first few cases are Ei = 3c/ L for M = 1 {N = 3) with Ei = 3c/ L 
and El = c(5 - VlO)/L, E2 = c(5 + VTO)/L for M = 2 (A^ = 5). 

For both the even and odd cases our results indicate that the Bethe roots approach 
the origin proportional to ^/c as c — for fixed A^ and L. We have checked this directly 
via numerical solution of the Bethe equations 

We now address the precise correspondence between the boson and BCS models. 
In view of the ground state energy £q = J2f=i ^1 boson model, we can map 

two bosons with opposite momenta ±kj onto one Cooper pair of bosons with energy 
Ej = 2k'j. In such a way, we can show that the ground state of the boson model with 
weakly repulsive delta interactions can be described by the BCS boson pairing model 
with strong pairing interaction. In particular, we need to consider the boson pairing 
model Hg 

n 

n = Y,^m + 2g J2 (17) 

^ i,i'=i 

based on the su{l, 1) algebra. Here = |aja£ + ^ denotes the particle number operator 
at level i and A\ = ha]a\ creates a boson pair with zero angular momenta at the level j. 
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The single-particle operators and a satisfy Bose commutation relations. Essentially, 
the single-particle state can be its own time reversal state for bosons at same levels. 
Therefore at every energy level there is a single pair with degeneracy f2 = 1 jTHj. For 
the boson pairing model, the pair energies are real. In the above equation, n is the 
number of unblocked energy levels and g is the pairing interaction. The energy of the 
pairing model ()17p is given by 

m 

e j=i 

where the pair energies Ej satisfy Richardson equations of the form JH] 

Here m denotes the total number of boson pairs and rf^ = y -|- ^ is the effective pair 
degeneracy of single-particle level £, where vi denotes the number of unpaired particles 
in level ^. Linking to the boson pairing model ^Hl; here the degeneracy Vti = 1 and 
z/£ = 0. 

In terms of the Cooper pair energies Ei participating in the scattering process, the 
equations (fT^ are 

1 2 _ n' 

~Yg^j^^E,-E, ~~eI' ^ ' 

for i = 1, . . . , r. Here n' = ^ + 2(m — r), where n is the number of unblocked energy 
levels, m is the number of Cooper pairs and r is the number of nonvanishing roots to 
order g. Direct comparison between Eqs. (jlOj) and ()2()|1 gives g = j^,r = m = M and 
n = 1, leading to a ground state of the BCS boson pairing model ()17|) with energy 

r 

^BCS = 5Z^i = 2^7r(n' + r-l). (21) 

i=l 

In the strong coupling limit, all energy levels can collapse into one multiply-degenerate 
level. It is thus reasonable that n = 1. Indeed, this state appears to lie in the physical 
regime of the boson pairing model ()17|1 . The ground state of the weakly interacting 
boson model thus corresponds to the ground state of the strongly coupled BCS boson 
pairing model. From equations (fTIH) and (fTH|) . the ground state energy per particle of 
the weakly interacting Bose gas follows as 

^0 _ V _ 2gBcs _ cjN - 1) 

This agrees with Eq. (jH)) and the results of jT] (see also Refs. j2Zll2H])- We see that it 
also coincides with the ground state energy per particle of the strongly coupled BCS 
boson pairing model ()17p . 

Further, if the Cooper pair is defined [29^ as 6kK = CfcaiCfc^i, b]^-j^ = c^jfC^.^! where 
k = ^{ki — K2) and K = (/ci + K2), it is shown that [&kK, ^u'r] = for k 7^ k', while 
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[^kK, ^Ic'k] = [^kK' ^I'k] ~ always hold. The single operators and c satisfy Fermi 
anticommutation relations. If the relative momentum vector of paired electrons lies 
inside the overlap of the two spherical shells in momentum space, the Cooper pair with 
nonzero centre of mass momentum, i.e. K = (fci + is possible and physical |29j . 
This means that the Cooper pairs can occupy a state of energy ex with same K but 
different values of k. Thus such Cooper pairs act as bosons and can collapse into a 
Bose condensate. In this BCS pairing model, the degeneracy of single-particle level fi^ 
would not be restricted to be even such that the quantum number n' in the equation 
flll|) can be fractional and negative. In this case, the single-particle states for a Cooper 
pair are not time-reversal eigenstates. For the standard BCS model, K = 0, such that 
the state k| is occupied so is the state — k|. Therefore we can expect that the BCS 
model of fermion pairs with nonzero centre-mass momentum would exhibit some of the 
behaviour of the ID boson gas with two-body delta interactions. 

In conclusion we have looked in detail at the many-body solution of the Bethe 
equations for the interacting boson model in the strongly and weakly repulsive coupling 
regimes. We found a remarkable connection with the Richardson equations for the 
strong coupling limit of the boson pairing model, suggesting a direct link between 
strongly interacting Cooper pairs and weakly repulsive bosons. This link deserves further 
investigation, particularly in light of the revival of interest in the interacting Bose gas 
and the BEC/BCS crossover 
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